Magnetization reversal of ferromagnetic nanodisk placed above a superconductor 
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Using numerical simulation we have studied a magnetization distribution and a process of mag- 
netization reversal in nanoscale magnets placed above a superconductor plane. In order to consider 
an influence of superconductor on magnetization distribution in the nanomagnet we have used Lon- 
don approximation. We have found that for usual values of London penetration depth the ground 
state magnetization is mostly unchanged. But at the same time the fields of vortex nucleation and 
annihilation change significantly: the interval where vortex is stable enlarges on 100-200 Oe for the 
particle above the superconductor. Such fields are experimentally observable so there is a possibility 
of some practical applications of this effect. 



PACS numbers: 75.75.+a, 75.60.Jk, 74.81.-g 



I. INTRODUCTION 

In the past few years considerable attention has been 
devoted to the investigations of magnetism in small nano- 
sized ferromagnetic particles. Such interest is caused by 
the opportunities for creating recording devices^ - — and 
ultrasmall magnetic field sensors^ based on the proper- 
ties of ferromagnetic particles. It is now well-understood 
that magnetization distribution in a single particle is de- 
termined by the competition between the magnetostatic 
and exchange energies. If a particle is small, it is uni- 
formly magnetized and if its size is large enough a non- 
uniform(vortex) magnetization is more energy prefer- 
able (see, for example, Refs. 0, 0, IE 0)- Besides the 
geometrical form and size, the state of the particle de- 
pends on many other factors. For example, by applying 
a homogeneous magnetic field we can cause nucleation or 
annihilation of the vortex. In an array of particles their 
magnetostatic interaction may have a strong influence on 
a particle magnetization. If the distance between parti- 
cles is rather small then magnetostatic interaction has a 
strong destabilizing effect on the vortex state, leading to 
a significant decrease in both the nucleation and annihi- 
lation fields^. 

The interplay between the ferromagnetism and super- 
conductivity can also lead to changes in the magnetiza- 
tion distribution. It was shown that in a ferromagnetic 
film put on a superconducting substrate the size of the 
domains is up to y/1.5 times smaller than for a film with- 
out a superconducting substrata^*— The experimen- 
tal investigations revealed changes in the magnetic field 
around Al/Ni submicron structures with a decrease in the 
temperature to values below T c . This effect was referred 
to the expulsion of the magnetic field by the supercon- 
ducting part of the Ferromagnetic/Superconductor (FS) 
hybrid structures— 

All investigations of the FS interaction deal with the 
changes of the ferromagnetic domain structure. There- 
fore the results thereof cannot be applied to the single- 
domain nanoparticles. On the one hand, the magneti- 
zation of a nanoparticle is more simple than the domain 



structure of a macroscopic ferromagnetic, therefore, the- 
oretical findings could be proved by experiments with 
nanoparticles. But on the other hand, the magnitude of 
the interaction between the nanoparticles and the super- 
conductor is strongly reduced due to the large values of 
the London penetration depth. 

In this work we investigate the phase transition be- 
tween the single-domain and the vortex state, and the 
process of magnetization reversal of a ferromagnetic 
nano-sized particle placed above the surface of a super- 
conductor. The superconducting state is described using 
the London approximation, that is, we assume that the 
particle cannot produce a vortex- antivortex pair. This 
assumption is not universally true, since if the particle 
dimensions are sufficiently large its magnetic field can de- 
stroy the Meijssner state of the superconductor—. The 
criteria of the London approximation applicability can 
be obtained in the following way. The largest magnetic 
field is generated by a single-domain particle. Then, the 
maximal magnetic flux through the surface of a super- 
conductor is 4> m = TtRhM , where R and h are the radius 
and the height of the particle and M is the magnetiza- 
tion. The flux <j) m must be less than the flux quantum 
4>o = 2, 07 • 10~ 7 gs ■ sm 2 in order to exclude the possi- 
bility of a vortex penetration. We consider a ferromag- 
netic particle with a saturation magnetization of about 
M s = 800 oe. Then we have the following limitation on 
the particle dimensions: Rh < 8- 10 3 nm 2 . So we consider 
only the particles which obey this condition. 

The equilibrium distribution of magnetization which 
gives a minimum to the energy functional is found by 
numerical simulation. For numerical simulation we use 
the approach based on the Landau-Lifshitz-Gilbert(LLG) 
equation for the dynamics of magnetic moments. This 
approach enables us to investigate metastable states 
which realize the local minimums of the energy func- 
tional. The metastable states are of great importance 
since the finiteness of the nucleation and annihilation 
fields values is the consequence of the energy barrier be- 
tween single-domain and vortex-like states. We obtain 
a phase diagram of the particle in the height/diameter 
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FIG. 1: Ferromagnetic nanoparticle above the surface of su- 
perconductor. 



plane for the transition from the single-domain to the 
vortex-like state. We find that for the realistic value 
of the London penetration depth the FS interaction is 
rather weak. The energy of the FS interaction is about 
100 times smaller than the self energy of the particle. 
Therefore, the FS interaction has no influence on the 
phase diagram. Nevertheless, the magnetization curve of 
the particle in an external homogeneous magnetic field 
changes significantly. Even the 1/100 energy shift means 
that the superconducting screening current produces a 
magnetic field of about 47rM/100 ~ lOOOe which leads 
to experimentally observable decrease of the vortex nu- 
cleation field and increase of the annihilation field. The 
paper is organized as follows. In Sec. 2 we derive the ana- 
lytical expression for the energy functional of our system. 
In Sec. 3 we show the most important features of the nu- 
merical simulation. In Sec. 4 we discuss the results and 
also propose possible experiments. Finally, the summary 
is given in Sec. 5. 



II. ENERGY FUNCTIONAL 



where J is the constant of exchange interaction. The 
second term E m in Eq. (JJJ is the magnetostatic energj*i&: 

E m = -^jn-Md 3 r. (3) 

Here H is the sum of magnetic field H m induced by the 
particle and magnetic field H s generated by the super- 
conducting current: 

H = H m + H s . 

The last term E ext in Eq.QJ is the energy of the inter- 
action between the magnetic moment and the external 
field H : 

E ext = - [ H Md 3 r. (4) 
Jv 

In order to obtain the energy as the functional of M we 
should express the magnetic field H in terms of M. To 
this end we should solve the Maxwell-London equations 
where the source is a magnetic current j m = cV x M. 

In the space outside the superconductor the system is 
described by the following equations: 

47T 

V x B = — j m (5) 



V ■ B = (6) 

where we introduce the vector of magnetic induction 
B = H + 47rM. The superconductor is described by 
the equations: 



V x B = — j, (7) 



The system we consider is a ferromagnetic particle 
placed above the surface of a superconductor (Fig^}. 
We assume that the superconductor occupies the whole 
half-space, so it has only one boundary plane, z = 0. 
The ferromagnetic particle is assumed to be made of soft 
magnetic material, so the energy of anisotrophy is equal 
to zero. The total energy of the system of particle and 
superconductor reads: 



E — E e + E m + E^t. (1) 

The first term E e is the energy of the exchange interac- 
tion: 

Ee = -2^2 / VM *| 2 + \^My\ 2 + \\7M Z \ 2 ) d 3 r, (2) 



B = -A£— Vxj s , (8) 

where j s is the superconducting current and Xi is the 
London penetration depth. For the boundary condi- 
tions we take the continuity of the magnetic induction 
B and its derivatives at the boundary plane z = 0. Since 
the Maxwell-London equations are linear, we first can 
consider a single magnetic dipole placed at the point 
fo = {xo,Uo,zo) above the surface of the superconduc- 
tor: 

M = mS(r - r ), 

where m is the magnetic dipole moment. Once this 
problem is solved, we can construct the solution for any 
magnetization distribution. After a bit troublesome but 
straightforward computation, we obtain the solution to 
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the single-dipole problem at the half-space z > in the 
form: 

B = -V0 + 4ttM, 
where the magnetic potential <j> is: 

4> = 4>d + 4>refl + 4>\- (9) 

The first term (f>d in Eq.© is the potential of a single 
magnetic dipole without a superconductor: 

cj) d (r) = -m- V, — - — 

The second term (j) re fi is the potential of the magnetic 
dipole reflected at the plane z = 0: 

<t>refi(r) = -m'-7 — — - , 
I 1 I ol 

where 

m* = (m x ,m y , -m z ). 

The third term <f>\ in Eq. © appears due to the finiteness 
of the London penetration depth value A^. The expres- 
sion for <px is: 



Here the first term is an ordinary dipole matrix: 

d 2 1 

Ddx 1 x 2 ( r ) = 7J a 1 fl 

ox\OXi |r| 

and the second term appears due to the presence of the 
superconductor 

where indexes x\ and x-i denote coordinates x, y, z and 
S z ,x 2 is the Kroneker's symbol. 

Let us now consider continuous distribution of the 
magnetic moment M(r). The magnetic field generated 
by the superconducting current is given as a sum of 
single-dipole solutions: 

H s = / D s (r-r*)M(r )d 3 r . 
Jv 

For the magnetic field H m generated by the magnetic 
current we cannot take the sum of the single-dipole fields 
since the dipole matrix D<j(r) has a nonintegrable singu- 
larity at point r = 0. Therefore, the expression for H m 
that we use is: 

H m (r)=-/ V-M(r ) r ~ r ° rf 3 r + 
Jv |r-r | 3 



A (r) = -m*-V.g(r-r*), (10) 

where 

g(r) = \ 2 L 2k 2 exp(-kz)J (kp)dk- 
Jo 

poo 

\\ \ 2k(k 2 + XI 2 ) 1 ' 2 exp{-kz)J (kp)dk, 
Jo 

p=(x 2 +y 2 y/ 2 . 

Here Jq(x) is the Bessel function. Then the single dipole 
magnetic field above the surface of the superconductor is 
expressed as: 

B(r) =H(r)+4^m,5(r-r ), 

H(r) = D(r, ro)m 
where D is the modified dipole matrix: 

D(r,r ) = D d (r-r )+D s (r-rS), (ff) 



/ (M(r ).n,(r )) . r \ d 3 r . 
Jdv \r-r \ 6 

Here, V, dV are the volume and surface of a particle and 
n s is the unit vector of the external normal to the surface 
at a current point. 

Finally, we consider the energy functional consisting 
of the self-energy of particle E , the energy of interac- 
tion with the superconductor Ei nt and with external field 

E e xt- 

E = E a + E int + E ext . 

Self-energy Eq includes exchange E e energy and the part 
Ed of magnetostatic energy E m : 

Eo = E e + Ed, 

where 

E d = -^J^U m -Md 3 r = 

1 f d 3 rd 3 r' 
-i / ( V .M(r))(V-M(r'))^r 
z JVxV |r — r i 
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1 f rf 3 rrf 3 r' 
+ i / ( V .M(r'))(M(r).n s )^L 

L JVxdV \ r — r I 

i r d 3 rd 3 r' 

— / (M(r')-n' s )(M(r).n s ) T — -. 

The energy of interaction £7j n t of the particle with the 
superconductor is the energy of the magnetic moment in 
the field produced by the superconducting current, and 
it is the other part of the magnetostatic energy E m : 

E ** = -\j H S 'Md 3 r = 

-\ [ M(r).D s (r,ro)M(ro)d 3 rd 3 r 

The energy E ext is the same as in Eq.Q and is given by 
Eq.Q. 

III. NUMERICAL SIMULATION 

To perform numerical calculation we use the same ap- 
proach as in Refs. Hand Hi The basis of our numerical 
simulation is the LLG equation for magnetization M(r,t) 
of a particle in the form 

^ = ^M,H eff l - [M, [M.H eff ll, 

(12) 

where 7 is the gyromagnetic ratio, a is the dimensionless 
damping parameter and t is time. The effective field 
H e / f is a variation derivative of the energy functional: 



Me// = ~SM 

The important feature of the LLG equation is that it 
describes the evolution of the magnetization distribution 
to the equilibrium. By varying the initial conditions we 
find different equilibrium states of our system. Then we 
choose the equilibrium state with the minimal energy to 
obtain the ground state. Choosing as the initial state the 
vortex-like or single-domain distribution and varying the 
external field H we determine the vortex annihilation or 
nucleation field as the critical field of transformation of 
the vortex-like to single-domain state or vice versa. 

To avoid a three-dimensional grid problem to be solved 
which needs large computer resources we assume that 
magnetization of a cylindrical particle does not depend 
on the coordinate z along the cylindrical axis. Then we 
integrate the relations for the energy over z and zq and 
obtain the energy as a functional of the magnetization 
which is a function of only two space variables. Then we 
define the effective field H e * f as a variation derivative of 



the obtained functional. The effective field does not de- 
pend on z either, so we have a three-dimensional problem 
reduced to the two-dimensional one. 

To develop a numerical method we divide the parti- 
cle into rectangular parallelepipeds with a square base 
of size a in the plane (x, y) and of height h and ob- 
tain approximate expressions for different parts of the 
energy functional using the grid values of magnetization 
M=(M x (p),M y (p),M z (p)). 

The expression for the magnetostatic energy E4 reads: 



E d [M] = -y E M (?') ' ^AP p')M(p) + E% (13) 

where p — (x,y) are the points of the square grid with 
step a on the (x, y) plane and matrix is the dipole 
matrix D f j integrated over the z coordinate: 

i> h d ( P )= [ dz [ z i) d (p,z r )dz'. 

JO J-z 

The additional term E® in Eq. 1)13(1 appears as the contri- 
bution of self-interaction in the cell. It depends only on 
the value of magnetization which is assumed constant, so 
it does not influence the effective field. 

The expression for the exchange energy reads: 

S e [M] = -^E^|M(p)-M(p')| 2 . (14) 

s p p' 

The internal summation in Ea. (|14[) is taken over all 
neighbors p' of the point p. 

The energy E ext and the interaction energy Ei„ t are, 
respectively: 

E elt ^a 2 ^H -M(p); (15) 
p 



E mt [M] = -y E M ^') ' 6J(p- p')M(p), (16) 
where matrix is 

rd-\-h rz+d+h 

i> h s (p) = dz i> s ( P ,z')dz'. 

Jd Jz+d 

We choose the cell size considering two factors. On 
one hand the size of the cell should be smaller than the 
characteristic exchange interaction length (a/ J/M 2 ) in 
order to describe the inhomogeneous magnetization cor- 
rectly. On the other hand, we cannot choose it very small 
because of the computation time limitations. 
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IV. RESULTS AND DISCUSSION 

For computation we choose the following parameters: 
the saturation magnetization M s = 800Oe, the exchange 
interaction constant J = 10~ 6 erg/ sm. The cell size is 
3nm x 3mn, while the exchange interaction length is ap- 
proximately 13rtm. The distance d and the London pen- 
etration depth have a strong influence on the interaction 
between a particle and a superconductor. For supercon- 
ductors with Xl ~ 50 — lOOnm the interaction should be 
largely reduced, since such values of Xl are comparable to 
the size of the particle. To compensate for this reduction 
we can make the distance d small. We choose A = 50nm 
and d = 5nm since such values are experimentally ob- 
tainable and make the effect of a particle-superconductor 
interaction quite distinctive. 

First of all, we investigate how the energy of the in- 
teraction between the particle and the superconductor 
depends on the value of Xl- We calculate the interaction 
energy Emt for the particle with a single-domain and 
vortex-like magnetization(Fig[2Jl. The dimensions of the 
particle are: h = lOnm and D = 50nm. The self-energy 
of this particle is Eq = 1.2122 • \Q~ ll erg for the single 
domain state and E = 1.8856 • 10~ 11 erg for the vortex 
one. The interaction energy decreases rapidly when Xl 
changes from to 30 — AOnm. For larger values of the 
interaction energy tends to zero asymptotically. At prac- 
tically interesting values Xl ~ 50 — lOOnm the interaction 
energy is about Ei nt ~ 0.0268 • 10~ 11 erg for the single- 
domain magnetization and Ei n t ~ 0.0014 • 10~ 11 erg for 
the vortex-like magnetization. Thus the interaction with 
a superconductor is much stronger for a single-domain 
particle. This is easily understood because the vortex-like 
magnetization produces a much weaker magnetic field 
than the single-domain magnetization. The interaction 
energy for the real values of Xl ~ 50nm is much smaller 
than the particle self energy: for the single-domain mag- 
netization it is Ei n t/Eo ~ 0.02. 

4 E int *10",erg 
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FIG. 2: The energy of interaction of a single-domain (filled 
circles) and vortex-like (open circles) magnetized particle with 
a superconductor as a function of Al- 

Then we study the ground state of the particle inter- 



acting with the superconductor. We define the ground 
state as a stable state with the lowest energy. There 
are two possible stable states for the particle: vortex-like 
and single-domain states. We obtain a phase diagram 
(see Fig|3J) where in the area above the phase boundary 
the ground state is vortex-like and below the boundary 
it is single-domain. Different curves on FigOU correspond 
to different values of Xl- 
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FIG. 3: Phase diagram of the vortex-single domain state 
transition. Filled circles- Xl = oo, open circles- Xl = 0, 
triangles- Xl — lOnm, rhombus- Xl — 50nm. 

Each phase boundary on Fig|3] has basic properties 
coming from the origin of phase transitions in a ferromag- 
netic nanoparticlei^. They are determined by the inter- 
play between the magnetostatic and exchange energy. In 
relatively large particles the ground state is vortex-like. 
When the dimensions (diameter D or height h) of the par- 
ticle are reduced then the single domain state becomes 
the ground one. The superconductor influence shows up 
through an enlargement of the area of small-sized parti- 
cles having the vortex-like state. It is clear that the shift 
of the phase boundary in comparison to the case A = oo 
depends on the ratio of the interaction energy to the self- 
energy of the single-domain particle: Ei n t/Eo. Since this 
ratio is very small Ei nt /Eo = 0.02 < 1 at A; = 50nm, 
the phase boundary remains almost unchanged. 

But the phase boundary between the two ground states 
is not important because at the large area around this 
curve on the phase diagram both the single domain and 
the vortex-like states are stable. When we cross the phase 
boundary from a large size to a smaller one, the vortex- 
like state becomes metastable, i.e., it is not the ground 
state but it does not transform to the single-domain state. 
This is also true for the single-domain state. The example 
of the area of metastability is shown on FigQ] 

Thus it is more important to investigate the bound- 
ary between the area of metastability and the absolute 
instability of the vortex-like or the single-domain state. 
This boundary depends on the external magnetic field 
Ho. It was found out that applying the in-plane field 
Ho we can cause nucleation of the vortex, i.e., make the 
single-domain state unstable, annihilation of the vortex, 
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FIG. 4: Example of the metastability area (between the 
curves) of the vortex-like state. The upper curve (filled cir- 
cles) is the boundary of the vortex-like state stability region, 
the lower curve(open circles) is the phase boundary between 
the vortex-like and single-domain ground states. 



i.e. make the vortex-like state unstable 8-9 . Practically, 
it is more convenient to vary the external magnetic field 
than the dimensions of the particle. That is why we 
investigate the vortex annihilation, H ann (FigJSJ) and nu- 
cleation, H nuc i (FigEJ) fields as functions of particle di- 
ameter D. The height of the particle is fixed: h = 20nm. 
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FIG. 5: The critical field H an n of vortex annihilation as a 
function of diameter D. The height of the particle h = 20nm. 
The curve marked with filled circles is for \l = oo, with open 
circles- for A = 50nm. 

The nucleation and annihilation of the vortex are parts 
of the process of magnetization reversal in a particle 
driven by the external magnetic field. If the particle is 
initially found in the single-domain state, then by ap- 
plying external field Ho in the opposite to the magnetic 
moment direction we first stimulate a transition to the 
vortex-like state, i.e., nucleation of the vortex. When we 
increase the external field further, the vortex annihilates 
and the particle comes to the single-domain state again 
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FIG. 6: The critical field H nuc i of vortex nucleation as a 
function of diameter D. The height of the particle h — 20nm. 
The curve marked with filled circles is for A^ = oo, with 
empty circles- for A = 50nra. 



but with a reversed magnetic moment. Note that the 
annihilation field H ann as always positive for the par- 
ticle with h = 20nm (see Fig|3J). The nucleation field 
Hnuci becomes negative when the diameter D is large 
enough(see FigEJ) • This means that the single domain 
state is unstable and the external field should be applied 
in order to prevent the nucleation of the vortex-like state. 
In more detail we investigate the magnetization reversal 
for a particle of D = lOOnm and height h — 20nm. At 
the zero field the particle is in the single-domain state 
with an average magnetic moment directed along the x 
axis. The external field is applied in the opposite direc- 
tion to the initial magnetic moment. The magnetization 
curve for this particle is shown on Fig[7|d). This curve 
describes the dependence of the average x-component of 
the magnetic moment M x on the external magnetic field. 
The distribution of magnetization at the stages of the 
magnetization reversal process is shown on Fig[7f2). 

The influence of a superconductor shows up in the re- 
duced value of the average magnetic moment, because it 
decreases the interaction energy Ei nt . Therefore, the nu- 
cleation field for the particle placed above the supercon- 
ductor is smaller and the annihilation field is larger than 
for the particle in the absence of superconductor. The 
difference between these fields increases with the size of 
the particle(see FigEl EH- When the diameter is quite 
large: D > lOOnm, the shift of the nucleation and anni- 
hilation fields is about 100 — 200oe. As we have noted 
before, this shift should be of order of the field produced 
by the superconducting current B s . To verify our re- 
sults we find < H s > averaged over the z coordinate: 
< H s > (x,y) = L H s (x,y, z)dz/h. For an example 
we take a single-domain particle with D = lOOnm and 
h = 20nm. On Fig|Hfa),(b),(c) we show the components 
of < H s > inside the particle. The magnetic field is nor- 
malized to the saturation magnetization M s — 800Oe. 

The x-component of magnetic field < H s > is the 
largest and it varies from 0.12Af s w lOOOe at the center 
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FIG. 7: (1) The process of magnetization reversal for Xl — oo 
(open circles) and Xl = 50nm (filled circles). The particle di- 
mensions are h = 20nm, D — lOOnm, (a),(d)-single-domain 
states, (b),(c)-vortex-like state. (2) The evolution of mag- 
netization distribution during the process of magnetization 
reversal. 




(b) 



(d) 



FIG. 8: (a)-(c) The level curves for the components < H 3X >, 
< H sy >, < Hsz > of the magnetic field produced by the 
superconducting current, (d)-magnetization of the particle. 
D = 150nm, h — 20nm, Xl — 50nm. 



ature T c we will have well observable shifts AH nuc i and 
AH ann of about 100 - 200Oe. Moreover, the shift of the 
annihilation field AH ann enables us to realize a selective 
magnetization reversal in an array of ferromagnetic par- 
ticles placed above the superconductor. If we neglect the 
interparticle interaction, then the magnetization curve of 
each particle is like the one shown on Fig0 

Let us assume that initially all particles in the array are 
in the vortex-like state. Applying the magnetic field of 
magnitude between H ann (T < T c ) and H ann (T > T c )(see 
Fig Hi)) will l eac l to a reversible displacement of the vor- 
tex core within each particle (see Fig|!j2a)). If then we 
destroy the superconductivity around one of the parti- 
cles, it will immediately fall into the single-domain state 
(see Fig(5f2b)). By removing the magnetic field we will 
have all particles get back to the initial state except for 
the one which will remain in the single-domain state. 
Thus we can operate with a single particle in an array 
without disturbing the state of other particles. 



of the particle to 0.06A/, » 50Oe at the edges. Thus 
AH nuc i and AH ann should be about 50 ~ lOOOe. Ac- 
cording to the results of simulation, the shift of the nucle- 
ation field for the cylindrical particle with D = 150nm is 
AH nuc i — 110Oe(see FigEl and the shift of the annihila- 
tion field is AH nuc i = — 120Oe(see FigEJ. Therefore, our 
estimation gives the right order of AH ann and AH nuc i. 

Let us now consider some possible experimental in- 
vestigations based on the effects which we have de- 
scribed. Placing a ferromagnetic particle with diame- 
ter D = 100 — 200nm and height h = 20nm above the 
superconductor and cooling it below its critical temper- 



V. CONCLUSION 

We presented the results of the numerical investigation 
of the magnetization reversal process in an external mag- 
netic field of a ferromagnetic particle placed above the 
surface of a superconductor. The numerical simulation 
is based on solving the Landau-Lifshitz-Gilbert equation 
for the dynamics of magnetic moment. The supercon- 
ductor is in the Meijssner state and the only parameter 
that affects the interaction between the particle and the 
superconductor is the London penetration depth A^. 

We have shown that for a realistic value of Xl = 50nm 
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the interaction energy is much smaller than the self- 
energy of the particle and the ground state of the parti- 
cles does not change significantly. But nevertheless, the 
magnetic field generated by the superconducting current 
leads to a decrease of the nucleation field and to an in- 
crease of the annihilation field by 100 — 200Oe. Based on 
the effect of the annihilation field shift when the super- 
conductor is cooled to temperatures below T c we describe 
the method of a selective magnetization reversal in an ar- 
ray of ferromagnetic particles. 
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FIG. 9: (l)-Part of the magnetization curve near the anni- 
hilation of the vortex. (2)-Selective magnetization reversal in 
an array of ferromagnetic particles. 
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